Full characterisation of polarisation states of light via direct measurement 
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Ascertaining the physical state of a system is vital in order to understand and predict its behaviour. 
Tomography is a standard approach used to determine the form of an unknown state. Here we show 
that an alternative approach, based on sequential weak and strong measurements, can be used to 
determine the density matrix in a simple, fast, and general style. We directly measure the probability 
amplitudes of a variety of pure polarisation states of light. We then generalise this experiment to 
directly measure the Dirac distribution and consequently determine the density matrix. Our work 
is the first to demonstrate the direct measurement of the full description of the state of a two- 
dimensional system, and it has applications to measurements in foundational quantum mechanics, 
quantum information, and physical chemistry. 



Introduction.— Measurement plays a vital role in the 
practice of science. This is especially so in the case 
of quantum mechanics, where the measurement process 
plays a fundamental role in the formulation of the the- 
ory. A crucial feature of quantum mechanics is that a 
measurement of a variable of a system erases informa- 
tion about the corresponding conjugate variable. The 
classic example is that determining the position of a par- 
ticle disturbs its momentum, and vice versa. These mea- 
surements, known as strong measurements, collapse the 
wavefunction such that no additional information can be 
obtained. 

In order to completely determine a quantum state, 
which is described in general by complex numbers, one 
must perform multiple measurements on many identical 
copies of the system. Quantum tomography [1 , proposed 
in 1957 [2], is one method of quantum state determina- 
tion that uses strong measurements [3H7]. Tomographic 
reconstruction entails estimating the complex numbers 
that describe the state from the real-valued probabili- 
ties that result from strong measurements. Consequently, 
this approach can be considered indirect state determi- 
nation due to the requirement of post-processing. 

The first demonstration of direct quantum wavefunc- 
tion measurement was recently reported [8]. In contrast 
to tomography, this method is considered direct because 
the measurement apparatus records the complex proba- 
bility amplitudes describing the state. A recent proposal 
outlined how to generalise the technique to measure all 
aspects of a quantum state [5]. Although familiar and 
convenient, the density matrix is not the only way to 
describe a general quantum state. A state can be ex- 
pressed in terms of its Dirac quasi-probability distribu- 
tion (or phase-space representative), which is informa- 
tionally equivalent to the density matrix p [TO] [11] . For 
the direct measurement of a quantum state, the Dirac 
distribution is a convenient representation as its elements 
can be written simply in terms of weak values. 

Directly measuring a quantum system relies on the 
technique of weak measurement: extracting so little in- 



formation from a single measurement that the state does 
not collapse [12-26 . The first measurement of a weak 
value was the amplified the transverse displacement be- 
tween the polarisation components of light induced by 
a birefringent crystal [15]. More recently, the technique 
was used to observe the transverse displacement a beam 
of light by only several angstroms [21] and an angular ro- 
tation on the order of femtoradians [22]. Weak measure- 
ment was recently proposed as a tool to study nonlinear 
optical phenomena with single photons by amplifying the 
apparent photon number [25] . Weak measurements have 
also allowed observation of apparent super-luminal veloc- 
ity [H] and the mapping of average photon trajectories 
after they pass through a double slit [24] . 

The main result of our paper is the implementation 
of a procedure for the direct measurement of the polar- 
isation state of coherent light. We demonstrate direct 
wavefunction measurement in a two-dimensional Hilbert 
space by weakly coupling the polarisation state of light to 
the spatial degree of freedom. We obtain the weak value 
by introducing a small spatial shift between the horizon- 
tal and vertical polarisation components, then strongly 
measuring the polarisation in the diagonal/anti-diagonal 
basis. Measurement of the mode's position and momen- 
tum is performed to read out the real and imaginary 
components, respectively, of the weak value. By means 
of the weak values, we report the direct measurement 
of the polarisation wavefunction [8]. In a separate ex- 
periment, by means of the weak values and associated 
probabilities, we report the direct measurement of the 
polarisation Dirac distribution [9]. We show that for the 
case of a pure state qubit, measuring a single weak value 
completely determines the state. The description of the 
polarisation state of light is the same in the classical and 
quantum regimes; thus, we choose to perform our exper- 
iment with classical light. This study extends previous 
work on polarisation weak measurements [I3j [15] [20] . 

Theory.— In any quantum measurement, the observer 
couples an unknown probe state to a pointer that reads 
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FIG. 1: Schematic representation of the experiment. State Preparation (a): The output of the single-mode fibre (SMF) is 
a near-Gaussian transverse mode of light. A polarising beamsplitter (PBS) and waveplate(s) create a known pure polarisation 
state. Weak Measurement (b): A quartz crystal at an oblique angle performs the weak measurement by introducing a 
small (compared to the beam waist) lateral displacement between horizontal and vertical polarisation components. Strong 
Measurement (c): A strong measurement in a basis mutually unbiased (diagonal/anti-diagonal) from the weak measurement 
is used to complete the direct measurement. Inset: To measure the wavefunction, a linear polariser oriented to transmit 
diagonally polarised light performs the strong measurement and post-selection; to measure the Dirac distribution, a A/2 wave- 
plate and calcite beam displacer perform the strong measurement. Readout (d): A 50:50 non-polarising beam splitter (NPBS) 
splits the light into two sub-ensembles. These are imaged in the near-field (dotted line, NF) and far-field (dashed line, FF) of 
the quartz crystal onto non-overlapping regions of the CCD camera. 



out the value. For example, a birefringent crystal can 
couple the polarisation state of light to the spatial degree 
of freedom; in this case, the initial polarisation state is 
called the probe state, and the spatial degree of freedom 
of the light is considered the pointer. If the pointer state 
is a Gaussian mode with width w, a strong measurement 
separates the eigenstates of the measurement operator 
A, by an amount S ^> w, such that the eigenstates are 
completely resolved. 



Weak measurements occur in the opposite regime, 
where the coupling is much less than the pointer width 
S <C w. In this case, the eigenstates of A are not resolved 
by the pointer, so the wavefunction does not collapse. 
Therefore, a subsequent measurement performed on the 
quantum state can be used to extract further informa- 
tion. If the subsequent measurement is strong, such that 
the eigenstates are resolved, we can choose to consider 
only the statistics of one particular outcome; this is called 
post-selection and the chosen outcome of interest is the 
post-selected state. The result of the weak measurement 



is called the weak value and is given by 
(A)J = 



. w _ {<p\Ap\4>) 
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where p is the density operator that describes the initial 
state and \<j>) is the final, post-selected, state [16J[T9]. In 
the case that the initial state is pure and may be de- 
scribed by the state vector (i.e., p = l^)^), the 
weak value in equation simplifies to the form first 
introduced by Aharonov, Albert, and Vaidman [T2] : 



(A) 



w 



(2) 



In the case that = |</>), the expectation value of the 
weak measurement is equal to the standard expectation 
value of the operator A. In general, the initial and final 
states may differ, and the weak value can be complex. 
For the specific case where the initial and final states are 
nearly orthogonal, the weak value can become arbitrarily 
large, leading to the amplification effect discussed above. 
The complex nature of the weak value, combined with 
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the fact that weak measurement does not significantly 
disturb the system, enables the direct measurement of 
the quantum state via weak measurements. 

The complex weak value is determined by character- 
ising the pointer. The pointer's position indicates the 
real part of the weak value Re[(A)^], and the pointer's 

momentum indicates the imaginary part Im[(A)^] [27] . 

In the specific case that the weak and final measure- 
ments are of mutually unbiased variables, the weak values 
have a direct relationship to the state description. The 
coefficients q of a wavevector \ip) that describes a pure 
quantum state can be written in terms of specific weak 
values: 

(k = (ai\ip) = y{^a x )f r (3) 

Here the weak measurement operator 7r a . = |a^)(a^| is 
the projector into the i th state of the basis A [8 . The 
factor v is a constant of normalisation independent of z, 
which may be taken to be real. Equation ([3| shows that 
the wavefunction describing a pure state can be directly 
measured by scanning weak measurements in basis A and 
post-selecting on a fixed state in the mutually unbiased 
basis 6, then normalising the wavefunction. 

An extension of the procedure that uses equation (|3| 
gives a technique to directly measure the most general 
description of the quantum state. The simplest such gen- 
eralisation entails measuring weakly in basis A, followed 
by recording the results of all outcomes of the strong mea- 
surement in basis B. In terms of the density operator p, 
the elements of the Dirac distribution that describes a 
general quantum state can be written in terms of specific 
weak values: 

S ij =(b j \a i )(a i \p\b j )=p bj (n ai )%. (4) 

That is to say, the (i, j) th element of the Dirac distribu- 
tion is equal to the result of the weak measurement of 
7r a . followed by post-selection on state bj, multiplied by 
the probability of successful post-selection p\>. = (bj\p\bj) 
[9]. Importantly, one can always invert equation Q and 
calculate the density matrix p from the measured Dirac 
distribution S. For further details on equations ([3| and 
Q, see the Supplementary Information. 

Experiment.— We perform two experiments. First, we 
implement the technique encapsulated by equation ([3| 
to measure a variety of pure polarisation wavefunctions. 
Second, we apply the technique encapsulated by equa- 
tion Q to measure the Dirac distribution of a variety of 
states. The only difference between the two experiments 
is in the nature of the strong measurement: in the first 
experiment, a single strong measurement outcome is re- 
quired; whereas in the second experiment, all eigenstates 
of the strong measurement are recorded. 

A brief summary of the experimental procedures is as 
follows; see Fig. [I] for a schematic. First, the probe (po- 
larisation) and pointer (spatial mode) states are prepared 



(a). Second, the weak measurement is performed by a 
quartz plate, which slightly displaces the two orthogonal 
polarisation components \H) and \V) of the probe lat- 
erally (b). Third, the strong measurement in the D/A 
basis is performed (c). To measure the wavefunction, we 
post-select the final state by projecting the polarisation 
into the diagonal state \D) using a linear polariser (LP) 
oriented to transmit diagonally polarised light. To mea- 
sure the Dirac distribution, a calcite crystal separates the 
components \D) and \A) so that they do not overlap. Fi- 
nally, the wavefunction or Dirac distribution is read out 
by imaging the near- and far-fields of the plane immedi- 
ately after the quartz onto separate regions of interest of 
a CCD camera (d). Two regions are used to read out the 
wavefunction and four are needed to read out the Dirac 
distribution. 

In order to demonstrate our ability to perform the di- 
rect measurement of the polarisation state, we measured 
the probability amplitudes of three sets of input polarisa- 
tion states, each corresponding to a different great circle 
on the Poincare sphere. The states were created by ap- 
propriate orientation of a half- wave plate and optionally, 
a quarter- wave plate. In the second experiment, we cre- 
ate a number of states in the same fashion and measure 
their Dirac distributions, then calculate the associated 
density matrix. 

Results.— Figure [2] shows the measured weak values and 
corresponding polarisation probability amplitudes as a 
function of input polarisation angle. Figure [3] shows the 
calculated Stokes parameters for each measured in 
Fig. [2] (blue points). We also show calculated Stokes pa- 
rameters for two additional paths around the Poincare 
sphere. Measured weak values and probability ampli- 
tudes are included in the supplementary materials for all 
these states. 

Figure [4] shows directly measured Dirac distributions 
and the corresponding density matrices of different po- 
larisation states. A variety of states are created for cali- 
bration as in the first experiment, but here all outcomes 
of the strong measurement are considered. 

Discussion.— Figure [2] shows that the largest divergence 
between theory and result in experiment 1 occurs when 
the initial state is anti-diagonal and therefore orthogo- 
nal to the post-selected state of diagonal. In this limit, 
the weak value is undefined, whereas the pointer reaches 
a maximum displacement (see [I2j [28]). This difficulty 
is overcome by the full state characterization technique 
performed in experiment 2. Recalling equation Q, it is 
precisely in this regime where p\ )j goes to zero, cancelling 
the effect of the breakdown of the weak-value approxi- 
mation. 

We note that density matrices determined by the tech- 
nique demonstrated herein are not guaranteed to be pre- 
cisely Hermit ian due to measurement noise. For example, 
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-67.5° (A) -45° (V) -22.5° (D) 0°(H) 22.5° (A) 
Half-wave plate angle (input polarisation state) 

FIG. 2: Results of experiment 1 with linearly polarised probe 
states, (a) Measured weak values plotted as a function of 
input polarisation angle. Error bars are shown only for the 
red points for clarity, (b) Real and imaginary components 
of the probability amplitudes determined by normalising the 
weak values of each test state, where = a\H) + f3\V). For 
both panels, the solid lines are the theoretical predictions of 
the real components, and the dotted lines are the theoretical 
predictions of the imaginary components. Inset: A Poincare 
sphere with the path taken indicated by the blue line. 




FIG. 3: Measured states on the Poincare sphere. Here we 
show a Poincare sphere with the set of directly measured 
states, indicated by their calculated Stokes parameters (<j x ) = 

(lp\7T D -7TA\lp), (cr y ) = (lp\7TR-7T L \lp}, (& z ) = (ip \tTH ~ \lp) • 

The blue points indicate states created by rotating the half- 
wave plate. The red (green) points indicate calculated Stokes 
parameters for states created by rotating the half- wave plate, 
followed by a quarter- wave plate at fixed angle 0° (45°). The 
solid lines indicate the paths taken for each data set. 
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FIG. 4: Results of experiment 2. The directly measured Dirac 
distributions and corresponding density matrices for the hor- 
izontal linear polarisation state \H) (a) and left-hand circular 
polarisation state \L) = 1/2 ((1 + i)\H) + (1 - i)\V)) (b). The 
axes of the Dirac distribution are the mutually unbiased bases 
H/V and D/A\ the density matrix describes the state in terms 
of only one basis H/V. 

the density matrices shown in Figure [4] have small imag- 
inary components along the diagonal, of magnitude on 
the order of the measurement uncertainty (« 3%). 

The similarity between equations ([3| and Q suggests 
a simple connection between the coefficients of the wave- 
function and the entries of its Dirac distribution. In the 
case that the state is pure, we may combine the two equa- 
tions to determine the real constant of proportionality 
that relates the two: 

c i Sij. (5) 

p bj 

We see that there is a column j of the Dirac distribu- 
tion that is proportional to the wavefunction q. It is, 
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in particular, the column corresponding to the choice of 
post-selection in equation Q that renders v independent 
of i (and hence can be taken to be real). 

Equation ([5| has particular relevance to our experi- 
ment for the states that have a constant of normalisation 
v equal to unity and a probability of post-selection equal 
to one half. These states lie on the great circle of the 
Poincare sphere that includes {\H), \R), |V), \L)} (red 
points in Fig. [3J weak values and probability amplitudes 
in Fig. [5] in the Supplementary information). Each state 
on this circle is from a basis that is mutually unbiased 
with respect to the strong measurement tt d . For these 
states, we see that q = 2Sij and hence the wavefunction 
is twice a column of the Dirac distribution. See Figure [4] 
for two examples. 

The technique we demonstrate compares favourably 
with quantum tomography with regards to inferring the 
density matrix from measurement results. Tomography 
via maximum likelihood estimation or least-squares fit- 
ting, which is an example of an inverse problem, becomes 
prohibitively difficult as the dimension of the state or 
number of particles in a multipartite state increases. The 
difficulty arises from the computational requirements of 
varying the vast number of fit parameters needed to esti- 
mate the state [7j [29] . In contrast, no fitting is required 
to determine the density matrix from the directly mea- 
sured Dirac distribution because it is calculated analyt- 
ically. Hence, we anticipate that for high-dimensional 
quantum systems especially, direct measurement will be- 
come a widely used technique for quantum state deter- 
mination. 

Conclusion.— In conclusion, we have performed direct 
measurements of general polarisation states of coherent 
light. We obtain our results through parallel measure- 
ments of the real and imaginary components of the weak 
value of polarisation. An important result is that a sin- 
gle weak value, corresponding to the weak measurement 
of only one observable, determines both complex coef- 
ficients of the pure state of a qubit. Our experiment 
was performed at classical light levels; however, this work 
paves the way for future direct measurement of discrete 
quantum systems. 

We anticipate that the general method prescribed in 
this work will be applied to a variety of two-level sys- 
tems. For example, the wavefunction or density matrix 
of polarisation-entangled photon pairs in quantum key 
distribution schemes and of spin qubits in quantum com- 
puting can be measured directly. If the calcite crystal in 
our second experiment is replaced with a second quartz 
plate as thin as the first, the density matrix of the ensem- 
ble can be determined with vanishingly small disturbance 
to its state before absorption by the camera. Indeed it 
should be possible, in principle, to measure the quantum 
state of an ensemble without destroying it, by weakly 
coupling the probe state of one particle to the pointer of 



another as in Ref. [20]. This has implications for quan- 
tum protocols, as well as physical chemistry, where one 
can imagine determining the quantum state of an ensem- 
ble without affecting it significantly. 

METHODS 

A near Gaussian pointer state is prepared by passing 
HeNe laser light through a single mode fiber (SMF). The 
probe is then prepared by polarising the light with a po- 
larising beam splitter (PBS) then rotating the angle of 
polarisation with a half- wave plate (A/2) and/or quarter- 
wave plate (A/4). 

The weak measurement is performed by coupling po- 
larisation information to the spatial degree of freedom 
of the light. Light incident at an oblique angle on an 
X-cut quartz crystal undergoes a polarisation-dependent 
parallel displacement. By aligning the extraordinary axis 
with the x-axis and rotating the crystal about the ?/-axis, 
horizontal and vertical polarisations become slightly sep- 
arated in x. We take the z-axis to be the direction of 
propagation of the light, x to be the transverse direction 
parallel to the optical table, and y to be the transverse 
direction perpendicular to the table. The angle of inci- 
dence was adjusted to ~ 40° to ensure the two optical 
paths through the crystal are equal, mod 27r, through 
the ordinary and extraordinary axes (based on crystal 
thickness of 700 jam). 

It is important that the pointer state be a Gaussian 
with a fiat phase- front. We collect the SMF output with 
a microscope objective (10 x), and focus the light onto the 
quartz crystal (~ 45 cm away). This ensures the phase- 
front is approximately flat over the region of interaction 
with the quartz. 

After strong measurement, the real part of the weak 
value is proportional to the average position (x) of the 
post-selected intensity distribution immediately behind 
the quartz. The quartz plane is imaged onto the camera 
by two sets of relay optics. The first set images (2/i- 
2/2 imaging system, fi = 100 mm, f 2 = 125 mm,) to 
a spatial filter (adjustable iris) that allows us to elimi- 
nate back-reflections created in the quartz crystal; the 
second set images (2/1-2/2 imaging system, /1 = 75 mm, 
fi = 250 mm,) the iris plane onto the camera. The imagi- 
nary part of the weak value is proportional to the average 
position of the intensity distribution in the far-field (p x ) 
of the quartz plane. A Fourier transform lens (/ = 300 
mm) maps the far-field distribution of the iris plane onto 
the camera. 

We establish the expectation value of each pointer by 
first integrating each intensity distribution I(x,y) along 
y to find I(x) = Y, y pixels I(x,y)Ay, followed by finding 
the average (x) = £^ pixels xI(x)Ax/ Ex pixels I(x)Ax. 
This procedure is repeated with the image of the far- 
field to establish (p x ), and for each strong measurement 
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outcome. 

The expectation values (x) and (p x ) of the pointer, and 
their corresponding standard deviations, are established 
by averaging 100 CCD images, each with a 2000-/is ex- 
posure time. The only exception is for the data used to 
calibrate the weak values for Figure [4}d, where we aver- 
aged over 50 CCD images, each with a 500-/is exposure 
time. This was to reduce the effect of spot drift over 
the course of the calibration run where many states are 
measured sequentially. 

A simple background subtraction is performed before 
calculating the pointer's position and momentum. We 
subtract the value of the minimum pixel from all pixels on 
each exposure, to reduce the effect that the background 
has on calculating the average. For the post-selection 
probability measurements used to determine the Dirac 
distribution, background subtraction is performed for 
each region of interest by subtracting the recorded inten- 
sity when the laser is blocked. The intensity after back- 
ground subtraction of the near-field image corresponding 
to the outcome \D) is Id and of \A) is I a- Thus the prob- 
abilities are calculated according to po = Id /(Id + I a) 
and p A = I a/ (Id + /a)- 

The weak value is obtained from average pixel number 

by 

(tt h )% = a(x)-b + i{c(p x )-d), (6) 

where a, 6, c, d are constants that must be determined 
by calibrating the measurement apparatus. Another set 
of calibration constants a'^b'^c'^d! must be determined 
for the post-selection of \A) to convert average pixel to 
(^hYa- We perform calibrations of the measurement ap- 
paratus by measuring the wavefunctions and Dirac dis- 
tributions of known pure states and comparing (x) and 
(p x ) to theoretically calculated weak values. 
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SUPPLEMENTARY INFORMATION 
Mutually unbiased bases 

Two bases A and B are said to be mutually unbi- 
ased if \(cti\bj}\ 2 = 1/N for all states \ai) in A and 
all states \bj) in S, where TV is the dimension of the 
Hilbert space in which lives and upon which p op- 
erates [30]. The indices i and j take the integers from 
to N — 1. For example, the horizontal/vertical polar- 
isation basis A = {\H}, \ V)} is mutually unbiased with 
respect to the diagonal/anti-diagonal polarisation basis 
B = {\D) = 1/V2(\H) + \V)), \A) = 1/V2QH) - \V))}. 



The wavefunction and weak values 

We summarize the technique proposed and performed 
by Lundeen et al in Ref. [8], where a spatial wavefunc- 
tion of an ensemble of single photons was measured di- 
rectly, using polarisation as the pointer. 

Recall that any state in the Hilbert space spanned 



by A can be written as 



(7) 



where each q = {a^iji) is a probability amplitude that 
cannot be accessed directly by strong measurements. The 
set {ci} is the wavefunction of the state expressed in the 
A basis. 

Consider a weak measurement of the projection opera- 
tor 7r a . = \<2i}(ai\, and post-selection on a state \<j>) = \bj), 
where \bj) and \ai) are states in mutually unbiased bases. 
In this case, equation gives the weak value 

(bj\ai)(ai\i 



<*0 



w 



v 



(8) 



where v is a constant of normalisation and is equal to 



(bj 



(9) 



Crucially, for the specific case where the weak and final 
measurements are mutually unbiased, the magnitude of 
the constant v is independent of i. In general, only for one 
choice of \bj) will the phase of v be independent of i. This 
ensures that the weak values are proportional through v 
to the complex probability amplitudes q. Thus, we can 
substitute the q in equation ^ with the normalised weak 
values from equation ^ to re-express the wavevector as 



(10) 



Observation of the wavefunction is therefore performed 
by measuring each coefficient in turn and normalising 
the magnitude of the wavefunction to unity, which deter- 
mines v. 

The wavefunction is determined by measuring the 
pointer state. The pointer's position indicates the real 
part of the weak value Re[(7r a .)^], and the pointer's mo- 



mentum indicates the imaginary part Im^Tr^)^ 1 



■ ./I), J EH- 

Returning to the example where the pointer is the spa- 



vin 



tial degree of freedom, the complex weak value (7r ai / b 
is established by measuring the expectation values (x) 
and (p x ) of the spatial mode: Re[(7r a .)^] = (x)/S and 
Im[(7r a .)^] = k(p x }/5. The constant of proportionality 
k = 2<j 2 /ft, with a the beam waist, is related to the mea- 
surement backaction [27 . 

In our work, we make use of the fact that not all the 
weak values (7r a .)^ are independent. Taking the sum of 
a complete set of weak values corresponding to projectors 
\ a i)( a i\ with fixed post-selection \bj), we obtain unity: 

N-l 



W 

'b 3 



{bj\ai)(ai\p\b 



i=0 




(bMbj) 



(ii) 
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This result is true for any post-selection \bj) = \<p) and 
relies on the fact that ^ |^)(a^| is equal to the iden- 
tity operator when the set of all |a$) is a complete basis. 
Therefore, measuring all but one weak value is sufficient 
to determine the wavefunction. Importantly, for a qubit, 
a single weak value determines the pure state in the two- 
dimensional Hilbert space. 

We now consider the direct measurement of the polar- 
isation wavefunction of a photon. Although the context 
of our work is polarisation, these results are general and 
apply to all two-state quantum systems. The polarisation 
state of a photon can be written as a linear combination 
of basis elements 



a\H)+p\V), 



(12) 



where a,/? are complex probability amplitudes, and 
\a\ 2 + \j3\ 2 = 1. Recalling equation ([3]), we may express 
the polarisation state as 

M=''((*H)%\H) + (ir Y )%\V)), (13) 

where the post-selected state is \D), and therefore mutu- 
ally unbiased with respect to the measurement operators 
tt h = \H)(H\ and tt v = \V)(V\. The constant v is sim- 
ply found by normalising the wavefunction and is equal 

to (KOST + \(k v )d\ 2 )~*- We now see that the wave - 

function is defined in terms of quantities that we can 
directly measure in the lab. Equation ( p~3| ) can be sim- 
plified by using the fact that (7r v )^ + (7r H )^ = 1 as in 
equation (11). Thus, we may eliminate (tt v )^, and the 
wavefunction is defined entirely in terms of (tt h )^ . 

The weak value (7r H )^ can be calculated from its defi- 
nition (equation 0) and the polarisation state (equation 



'D 



(D\H)(H[ 



a + /3 



cos 6 + e ltp sin 9 



(14) 



We take the convention here that the horizontal com- 
ponent has a real coefficient, such that the coefficients 
are defined by a = cos# and j3 = e^sin#, where 9 is 
the polarisation angle and ip is the phase difference be- 



tween components. The result of equation (14) can also 



be obtained through the methods of classical optics. The 
locations of the peaks of the intensity distributions ob- 
tained from the interference of two overlapping Gaussian 
modes, in the near- and far-fields, are proportional to 



the real and imaginary parts of equation (14), respec- 
tively [15]. In this representation, the two degrees of 
freedom of a polarisation qubit are 9 and p; in terms of 
weak values, they are the independent real and imagi- 
nary parts of (7r H )3f. Setting a convention such as this 
for the phase fixes the gauge, or reference frame, with re- 
spect to which the wavefunction is measured. Choosing 



v real actually sets a different phase convention, where 
the wavevector has an additional global phase coefficient 
exp(— iarg[a + /?]). We plot the wavefunctions in Fig. [5] 
using this phase convention, since it is the natural choice 
that comes from experiment. The theory curves of Fig- 
ures [5J[5j and [6] are calculated according to equation (14) 
by setting ip = 0,7r/2,7r/4 respectively, then setting the 
wavefunct ion's phase as discussed above. 



The Dirac distribution and weak values 

We summarize the technique proposed by Lundeen and 
Bamber in Ref. [9]. The Dirac distribution, named for 
Dirac who proposed it in Ref. [10 , is investigated thor- 
oughly by Chaturvedi et al in Ref. [11]. In our work, 
we use Chaturvedi et aZ.'s "left" phase-space represen- 
tative (Dirac distribution) throughout, and discuss only 
the discrete (i.e., N- level) Hilbert space version. 

The entry in the i th row and the j th column of the 
Dirac distribution of the state described by p is given by 



Sij = (bj\ai)(ai\p\bj). 



(15) 



Recalling the weak value defined in equation ([!]), for weak 
measurement of a projector ir a . , we obtain 



(16) 



Substituting equation ([16]) into equation ([15]), we now 
see that the elements of the Dirac distribution that de- 
scribes a general quantum state can be written in terms 
of particular weak values that can be directly measured: 



(17) 



That is to say, the (i,j) th element of the Dirac distribu- 
tion is equal to the result of the weak measurement of 
7r a . with a post-selection on state 6j, multiplied by the 
probability of successful post-selection = (bj\p\bj). 
The generalised technique uses the outcome of the strong 
measurement, and since a mixed state density matrix 
(and hence Dirac distribution) has more degrees of free- 
dom (N 2 — 1) than a pure state vector (27V — 2), the 
strong measurement projection must be scanned through 
the basis S, in addition to scanning the weak measure- 
ment through the complimentary basis A. 

To calculate the density matrix p from the Dirac distri- 
bution 5, a simple row- by-row discrete Fourier transform 
can be used, since the H/V basis and D/A basis are not 
only mutually unbiased but Fourier bases. The relation 
is 



N-l , 2 
k=0 



k)j 



(18) 
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In equation (18), we use i to denote the imaginary unit 



which is distinct from the index i. 

In terms of experiment outcomes, the Dirac distribu- 
tion is 

5= > Left basis axis (H\, (V\. 

\PoMn PaMa J J 
Right basis axis \D), \A) 

(19) 

Using equation (18) with N = 2, we obtain the density 
matrix 



FIG. 5: Results of experiment 1 with elliptical probe states 
(quarter- wave plate fast axis 0°). (a) Measured weak val- 
ues plotted as a function of half- wave plate angle. Error bars 
are shown only for the red points for clarity, (b) Real and 
imaginary components of the probability amplitudes deter- 
mined by normalising the weak values of each test state, where 
— a\H) + P\V). For both panels, the solid lines are the 
theoretical predictions of the real components, and the dotted 
lines are the theoretical predictions of the imaginary compo- 
nents. Inset: A Poincare sphere with the set of intended 
input states to directly measured indicated by the red line. 



P = 

\Pd(^v)d -Pa(^v)a Pd(^v)d +Pa{^v)j 

Both axes of this density matrix use the H/V basis. 
Additional Data 



(20) 
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Half-wave plate angle (input polarisation state) 
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FIG. 6: Results of experiment 1 with elliptical probe states 
(quarter- wave plate fast axis 45°). (a) Measured weak val- 
ues plotted as a function of half- wave plate angle. Error bars 
are shown only for the red points for clarity, (b) Real and 
imaginary components of the probability amplitudes deter- 
mined by normalising the weak values of each test state, where 
= a\H) + P\V). For both panels, the solid lines are the 
theoretical predictions of the real components, and the dotted 
lines are the theoretical predictions of the imaginary compo- 
nents. Inset: A Poincare sphere with the set of intended 
input states to directly measured indicated by the green line. 



